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ABSTRACT 

Context. Apart from being an important coolant, water is known to be a tracer of high- velocity molecular gas. Recent models predict 
relatively high abundances behind interstellar Shockwaves. The dynamical and physical conditions of the water emitting gas, however, 
are not fully understood yet. Using the Herschel Space Observatory, it is now possible to observe water emission from supersonic 
molecular outflows at high spectral and spatial resolution. Several molecular outflows from young stars are currently being observed 
as part of the WISH (Water In Star-forming regions with Herschel) key program. 

Aims. We aim to determine the abundance and distribution of water, its kinematics, and the physical conditions of the gas responsible 

for the water emission. The observed line profile shapes help us understand the dynamics in molecular outflows. 

Methods. We mapped the VLA 1623 outflow, in the ground-state transitions of o-H 2 0, with the HIFI and PACS instruments. We 

also present observations of higher energy transitions of o-H 2 and ;>H 2 obtained with HIFI and PACS towards selected outflow 

positions. From comparison with non-LTE radiative transfer calculations, we estimate the physical parameters of the water emitting 

regions. 

Results. The observed water emission line profiles vary over the mapped area. Spectral features and components, tracing gas in 
different excitation conditions, allow us to constrain the density and temperature of the gas. The water emission originates in a 
region where temperatures are comparable to that of the warm H 2 gas (Ti 200 K). Thus, the water emission traces a gas component 
significantly warmer than the gas responsible for the low-7 CO emission. The water column densities at the CO peak positions are 
low, i.e. N(H 2 0) = (0.03 - 10) x 10 14 cirT 2 . 

Conclusions. The water abundance with respect to H 2 in the extended outflow is estimated at .X(H 2 0) < 1 x 10~ 6 , significantly lower 
than what would be expected from most recent shock models. The H 2 emission traces a gas component moving at relatively high 
velocity compared to the low-/ CO emitting gas. However, other dynamical quantities such as the momentum rate, energy, and 
mechanical luminosity are estimated to be the same, independent of the molecular tracer used, CO or H 2 0. 

Key words. ISM: individual objects: VLA 1623 - ISM: molecules - ISM: abundances - ISM: jets and outflows - stars: formation - 
stars: winds, outflows 



1. Introduction 

Outflows are spectacular signposts of the star formation process, 
where swept up material from the parental cloud, excited in the 
interaction with the ambient medium, can be observed on par- 
sec scale distances from the central source (see e.g. Richer et al. 
2000; Tafalla & Bachiller 2011, for reviews). Outflows from 
young stars are most likely a consequence of a re-distribution 
of angular momentum during the process of proto-stellar accre- 
tion. Jets launched from a region close to the proto-stellar disk 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



(Banerjee & Pudritz 2006) likely drives them and shock waves 
interacting with the ambient medium decelerate the gas in the 
aftermath of the shock. Emission from H2O is produced in the 
post-shock medium where the molecules can cool on a relatively 
long time scale but the emission may also stem from within the 
jet itself (see e.g. Eisloffel et al. 2000). Momentum conservation 
implies that the bulk of the flow is moving at a moderate ve- 
locity and radial velocities should increase along the flow axis 
(see e.g. Cabrit et al. 1997). Previous spectrally resolved obser- 
vations reveal complicated line profiles where gas moving at dif- 
ferent velocities (with respect to the observer) show up as spec- 
tral features. The shape of the line profiles is a consequence of 
the geometry, abundance, and velocity fields but it also reflects 
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the optical depth in the observed lines. In some of the observed 
sources, high velocity components, hereafter referred to as bul- 
lets, show up (see e.g. Kristensen et al. 201 1). 

The water abundance in the outflow is expected to reflect 
the type of shock physics in action (see e.g. Bergin et al. 1998; 
Kaufman & Neufeld 1996; Hollenbach et al. 1989). Depending 
on the physical conditions in the flow, water abundances with 
respect to H2 can vary by several orders of magnitude in the 
post-shock gas, and can be as high as 10~ 4 (Flower & Pineau 
des Forets 2010). H2O line profiles and abundances are often 
compared to those of CO and other shock tracers, readily ob- 
served with ground-based telescopes. However, recent observa- 
tions with HIFI indicate that emission from low-/, CO and H2O 
does not always probe gas with comparable excitation condi- 
tions. Likewise, emission from these two molecules does not 
necessarily correlate spatially. Whilst the CO emission plausibly 
traces entrained gas, it has been found that H2 can be a valuable 
probe, tracing shocked gas at excitation conditions similar to that 
of H 2 (Santangelo et al. 2012; Nisini et al. 2010). 

The radiation from most water transitions is effectively 
blocked out by the Earth's atmosphere so these transitions can- 
not be observed using ground-based facilities. Instead, space 
based observatories such as SWAS (Melnick et al. 2000) and 
Odin (Nordh et al. 2003) have been used to observe the lowest 
rotational transitions of water. Using these facilities, the abun- 
dance of this molecule has been constrained (Franklin et al. 
2008; Bjerkeli et al. 2009). It has also been shown that the 
H2O abundance is generally lower than what was previously ex- 
pected. 

In this paper we present a study of the outflow powered 
by VLA 1623 carried out with the Herschel Space Observatory 
(Pilbratt et al. 2010). VLA 1623 is located in the constellation of 
Ophiuchus at a distance of 120 pc (e.g. Lombardi et al. 2008) 
near the p Oph A cloud core (Loren et al. 1990). This core is 
part of the star-forming cloud L1688, where a large number of 
stars have recently been formed (see e.g. Bontemps et al. 2001; 
Duchene et al. 2004). Apart from being located in one of the 
most nearby star-forming regions it is also considered to be the 
prototype of a very young stellar object, i.e. a Class source 
(Andre et al. 1990, 1993). The central source drives a collimated 
bipolar CO outflow extending more than 0.5 pc in the south- 
eastern direction (Andre et al. 1990, 1993; Dent et al. 1995). 
The nature of the source itself has in recent years been a matter 
of debate. For instance, Looney et al. (2000) suggested that it is 
a binary system with two circumstellar disks. This has, however, 
been challenged by the interferometric observations presented 
in Ward-Thompson et al. (2011). An alternative interpretation 
is that the observed continuum sources close to VLA 1623 are 
knots in the north-western jet (Maury & Andre 2012). 

Here we present HIFI (Heterodyne Instrument for the 
Far Infrared) and PACS (Photodetector Array Camera and 
Spectrometer) observations of the ground state transitions of or- 
f/10-water at 557 and 1670 GHz. Two maps were obtained, cov- 
ering a 5' x 2' region of the VLA 1623 outflow and the southern 
part of the p Oph A cloud core. Also, a line survey was con- 
ducted towards two outflow positions. Combined, these datasets 
provide us with kinematic and spatial information at high reso- 
lution. This allows us to constrain the physical and dynamical 
parameters of the gas responsible for the observed H2O emis- 
sion. 

The paper is organised as follows; the HIFI and PACS obser- 
vations are presented in Section 2 and the main results are pre- 
sented in Section 3. The obtained maps and line survey are dis- 
cussed in detail in Section 4 where we also perform an excitation 



analysis towards selected positions in the outflow. In Secction 5, 
our main conclusions are briefly summarised. 

2. Observations and data reduction 

The mapping observations of VLA 1623 were obtained with the 
HIFI (de Graauw et al. 2010) and PACS (Poglitsch et al. 2010) 
instruments aboard Herschel between 17 February 2010 and 12 
March 201 1. In addition, selected transitions of ortho- and para- 
water were observed towards two outflow positions. These ob- 
servations were carried out as part of the Water In Star-forming 
regions with Herschel (WISH) key program (van Dishoeck et al. 
2011). 

The data were processed using the Herschel Interactive 
Processing Environment (HIPE, Ott 2010). Baseline subtrac- 
tion and subsequent data reduction were performed in XS 1 and 
Matlab. 

2.1. HIFI 

2.1.1. HIFI - map 

A map of the region surrounding VLA 1623 was obtained in 
the 557 GHz ground state rotational transition of 0-H2O on 
1 September, 2010 (see Fig. 1). The 3.5 meter Cassegrain tele- 
scope has a Half Power Beam Width (HPBW) of 38" at 557 GHz 
and a pointing accuracy of the order of 1" (Garcfa-Lario et al. 

2011) . The calibration uncertainty is ~10% (Roelfsema et al. 

2012) and the main beam efficiency at 557 GHz has been mea- 
sured (from observations of Mars) to 0.76 (Olberg 2010). The 
data were calibrated using the Herschel Interactive Processing 
Environment (HIPE) version 7.2.0. In this paper, data from 
the Wide Band Spectrometer (WBS) are used. The WBS is 
an acousto-optical spectrometer with a bandwidth coverage of 
4 GHz. The spacing between the channels is 500 kHz which 
corresponds to 0.3 kms~' at 557 GHz. The HIFI map was ob- 
tained using the "On-The-Fly (OTF) Maps with Position-Switch 
Reference" observing mode, where the telescope is scanning the 
area. The reference spectrum was obtained toward a position 
which was located 14' away. In Fig. 1, however, the data are pre- 
sented on a rectangular grid. A Gaussian regridding procedure 
was used to compute the spectra at each position. The readout 
positions for the OTF map, used as input for the regridding, are 
presented in Fig. 2. Here, the positions for the horizontal (H) and 
vertical (V) positions were calculated separately. 

The Herschel HPBW at 557 GHz is about a factor of 3 
larger than that at 1670 GHz. To improve the spatial resolution, 
a Statistical Image Deconvolution (SID) method was used to de- 
convolve the HIFI map. This method is based on Bayes' theorem 
and is described in detail in Rydbeck (2008). The spatial resolu- 
tion of the HIFI map, after de-convolution, is estimated at ~20" 
in the region where the signal-to-noise is high. 

2.1.2. HIFI - pointed observations 

Selected transitions (see Table 1) were observed towards two 
outflow positions showing enhanced emisson in CO (2-1) (see 
Andre et al. 1990). These positions are referred to as Blue 1 
(Bl) and Red 2 (R2) hereafter and the offset with respect to 
VLA 1623 is (22", -11") for the Bl position and (-57", +29") 

1 data reduction software developed by P. Bergman at the 
Onsala Space Observatory, Sweden; http://www.chalmers.se/rss/oso- 
en/observations/ 
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Fig. 1. H 2 0(lio - loi) map obtained with HIFI. The data have been regridded onto a regular grid in RA and Dec. The spacing between each 
position is 19" and the map is centred on the Class source VLA 1623 at 112000 = 16 h 26 m 26.4!, foooo = -24°24'31". All spectra are plotted on the 
r m b-scale and velocities are with respect to the l>lsr (from -50 to +50 km s~'). Note that the spectra in the enlarged insets have not been baseline 
subtracted, but corrected such that the offset shows the single side-band continuum level. The positions of B 1 and R2 are indicated with coloured 
dots. 
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Fig. 2. Upper left panel: Colour map shows the integrated intensity over the H 2 0(lio - lot ) line from -40 to +55 kms -1 . Contours show the 
level of the 557 GHz continuum in steps of 14 mK where the first contour is at 79 mK. The HPBW is 38". Upper right panel: The colour map 
shows the flux of the H2O (2n - lot ) line. Contours show the continuum level in steps of 4 x 10 -4 erg cm" 2 s -1 yum -1 sr -1 where the first contour is 
at 4 x 10 -4 ergcirT 2 s -1 /jirT 1 sr -1 . The HPBW is 13". Lower left panel: The blue and red contours show the integrated H2O (lio - lot) intensities 
in the blue and red wings for the de-convolved map. Red contours are in steps of 0.85 K kms -1 and the first contour is at 0.16 K kms -1 . Blue 
contours are in steps of 1.14 K kms -1 and the first contour is at 0.16 K kms -1 . The readout positions for the OTF map are indicated with black 
dots. The horizontal and vertical polarisations are separated by 6" in the north-east to south-west direction. Lower right panel: The colour map 
shows the integrated intensity for the H 2 (2 l2 - lot) map convolved to 20" resolution. 

In all panels, the peak positions of the CO(2-l) emission are indicated with blue (Bl) and red (R2) dots. The star S 1, the young stellar object 
GSS 30 and other HH objects and near-infrared knots are indicated with labels and yellow squares. The beam sizes (or equivalent resolution for 
modified maps) are indicated with black circles in the lower right corners of each figure. 



for the R2 position. The observations were carried out in "Dual 
Beam Switch (DBS)" mode where the beam is switched by the 
internal chopper to a sky position 3' away from the target. 

2.2. PACS 

2.2.1. PACS - map 

A PACS map of the 1670 GHz ground state transition of 0-H2O 
was obtained on 1 September, 2010 (see Fig. 2). The map was 
obtained in line spectroscopy mode and the wavelength range 
was centred on the H2O (2i2 - loi) transition at 179 /on. A map 
of the VLA 1623 outflow was obtained by combining three raster 
maps centred at three different positions along the axis of the 
outflow. Each map consists of 3x3 PACS frames with a spatial 



separation of 40". The PACS instrument is a 5x5 integral field 
unit array having a spatial separation of 9'.'4 between each pixel. 
TheH20(2i2 - loi)lineis not resolved at R = 1700 correspond- 
ing to a FWHM of the line of ~175 kms -1 . The data were re- 
duced using Hipe v.6.0. The observed fluxes were normalised 
to the telescope background and subsequently converted to ab- 
solute fluxes using Neptune as a calibrator. The uncertainty of 
the flux calibration is ~30%. Software such as IDL and XS have 
been used to obtain continuum subtracted line maps. 

2.2.2. PACS - pointed observations 

Two ortho- and two para-transitions were observed using the 
PACS array in line spectroscopy mode (see Table 1) . These 
observations were carried out towards the blue- and red-shifted 
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Table 1. Observations carried out with PACS and HIFI. 



Instrument - Band 


Transition 




Spin 


Frequency 
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HPBW/FOV" 


Date 


f int 


Date 


? int 








isomer 


(GHz) 


(K) 


(") 


(YYMMDD) 


(sec) 


(YYMMDD) 


(sec) 


Mapping observations 






















HIFI-1 


H 2 O(l 10 - 


loi) 


ortho 


556.936 


61.0 


38 


100901 


6552 






PACS 


H 2 0(2 12 - 


lot) 


ortho 


1669.905 


114.4 


13 


100227 


11154 






oingl€ pointing 














R2 (-57", +29") 


Rl ( j-ll" 
Dl \-tLL ,- 


1 1 "\ 


HIFI-1 


H 2 O(l 10 - 


lot) 


ortho 


556.936 


61.0 


38 


100929 


878 


100901 


440 


HIFI-2 


H 2 0(2„- 


202) 


para 


752.033 


136.9 


28 


110312 


430 


100916 


454 


HIFI-4 


H 2 0(3 P - 


3o 3 ) 


ortho 


1097.365 


249.4 


19 


110217 


1658 


110217 


829 


HIFI-4 


H 2 0(1„- 


Ooo) 


para 


1113.343 


53.4 


19 


110217 


765 


110217 


829 


HIFI-6 


H 2 (2 12 - 


loi) 


ortho 


1669.905 


114.4 


13 


110218 


1008 


110218 


1058 


PACS 


H 2 O(3 03 - 


2,2) 


ortho 


1716.770 


196.8 


47 x47 


100227 


3098 


100226 


3098 


PACS 


H 2 0(3i3- 


202) 


para 


2164.132 


204.7 


47 x47 


100227 


3098 


100226 


3098 


PACS 


H 2 O(4 04 - 


3n) 


para 


2391.573 


319.5 


47 x47 


100227 


3098 


100226 


3098 


PACS 


H 2 Q(2 21 - 


lio) 


ortho 


2773.977 


194.1 


47 x47 


100227 


3098 


100226 


3098 



Notes. (a) For the PACS single pointing observations, the Field Of View (FOV) of the array is indicated in the table. 



outflow lobes with the central spaxel centred on the B 1 and R2 
positions respectively. For the excitation analysis presented in 
Sec. 4. 1 , line fluxes have been extracted in a circular region (32") 
centred on Bl and R2. Using this method, line fluxes should be 
accurate to within 30%. 



3. Results 

The H2(D(lio - lot ) and H20(2i2 - lot) lines are detected to- 
wards the outflow emanating from VLA 1623 (Fig. 1 & 2). In 
addition, several other transitions of 0-H2O and P-H2O are de- 
tected towards the CO peak positions, R2 and B 1 . We note, how- 
ever, that these positions do not correspond to peaks in the H2O 
emission. A summary of the observed H2O lines is presented in 
Table 2. 



3.1. HIFI map 

The line profile shapes observed towards individual positions 
change significantly throughout the H2OU10- lot) map ob- 
tained with HIFI. In Fig. 1, we present the profiles towards 114 
positions with five positions of particular interest enlarged. In 
summary, the line profiles exhibit 3 different components possi- 
bly tracing gas that originate in different physical regions. 

3.1 .1 . Absorption feature 

An absorption feature, which in some positions reaches be- 
low the observed continuum level, is visible across the map 
(^lsr - 3.7 kms~'). The absorption is detected at the majority 
of the positions in the map and is most prominent in a region 
centred (+25", +25") offset from the position of VLA 1623, i.e. 
near SMI (Motte et al. 1998). No clear, large-scale velocity dif- 
ferences for the absorption feature can be seen (see Fig. 3). The 
southern part of the core, where no or little outflowing gas is de- 
tected, is presented in spectrum 4 of Fig. 1 and shows hardly any 
absorption. 



^ 0.3 



1 1 1 1 1 




VLA1623 . 

GSS30 

R2 

Bl 










L 


. . 1 1 1 


1 


1 1 . . 



-40 -20 _, 20 40 



Fig. 3. H 2 0(lio - Iol ) spectra towards four selected positions in the 
map. No large-scale velocity gradients can be detected from the cen- 
troid velocity of the absorption feature. The absorption feature reaches 
below the continuum level where the continuum is strong, i.e. towards 
VLA 1623 and Bl. Note that the spectra not have been baseline sub- 
tracted. 

3.1.2. Spectral wings 

Towards the outflow, gas motions are clearly traced in the form 
of broad spectral wings. In the central region, both red-shifted 
and blue-shifted emission is detected. The blue-shifted part of 
the flow is free from contamination from other flows. Here, no 
red-shifted components are detected (see Fig. 1, spectrum 2). 
On the other hand, a blue-shifted component is clearly visible 
in the same region as the red part of the flow emanating from 
VLA 1623 (see Fig. 2 & spectrum 3 of Fig. 1). Here, the H 2 
emission line profiles show a strong blue-red asymmetry. The 
blue-shifted component is most likely associated with GSS30 
(see e.g. Tamura et al. 1990) and the peak of this emission is lo- 
cated close to the position of the Herbig-Haro object HH313 A 
(Davis & Eisloffel 1995). However, the presence of blue-shifted 
gas in the region close to R2 is not easily associated with GSS30. 
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Table 2. Integrated intensities of the lines observed with PACS and HIFI during the line survey. Numbers in parentheses are 1 a uncertainties. 











IX _ 




R 1 
1 




Transition 




Frequency 


Avlsr 


J T mb dv 


■* mb,rms 


/ ?mb dv 


T mb,rms 






(GHz) 


(km s" 1 ) 


(Kkms- 1 ) 


(mK) 


(Kkms- 1 ) 


(mK) 


H 2 O(l 10 - 


lot) 


556.936 


-30 to 30 


3.4 (0.03) 


8 


3.6 (0.05) 


12 


H 2 0(2„ - 


2 02 ) 


752.033 


-30 to 30 


0.3 (0.08) 


22 


0.8 (0.08) 


21 


H 2 0(3 12 - 


3 03 ) 


1097.365 






28 




30 


H 2 0(1„ - 


Ooo) 


1113.343 


-30 to 30 


1.4 (0.10) 


36 


1.9 (0.10) 


36 


H 2 0(2 12 - 


lot) 


1669.905 






106 




106 


H 2 O(3 03 - 


2 I2 ) 


1716.770 


-200 to 200 


0.7 


0.60 


0.6 


1.05 


H 2 0(3 13 - 


2 02 ) 


2164.132 


-200 to 200 


0.2 


0.19 




0.31 


H 2 O(4 04 - 


3o) 


2391.573 






0.19 




0.22 


H 2 0(2 21 - 


lio) 


2773.977 






0.29 




0.15 



Notes. <D) For the HIFI observations, the frequency bin size when calculating the rms, is the same as the channel spacing. For the PACS data, 
the frequency bin sizes at 1717, 2164, 2392, and 2774 GHz are 0.3, 0.5, 0.6, and 0.7 GHz respectively. The rms is in these cases measured in the 
central spaxel. 



In fact, blue-shifted emission is detected at all positions along 
the red part of the outflow. This could indicate that the north- 
western flow is nearly in the plane of the sky (Andre et al. 1990). 
At l<lsr < +20kms~ 1 , there is a trend where the velocity is in- 
creasing with distance from the central source (see Fig. 4). 



3.1.3. High-velocity gas 

Spectral features ("bullets") are detected towards several posi- 
tions in the map (see e.g. spectrum 2 in Fig 1). These features 
trace gas at velocities higher than those traced by the outflow 
wings. There is, however, no tendency for the highest velocity 
gas to be observed at a certain distance from the central source. A 
channel map of the region reveals high velocity gas at positions 
far from as well as close to the driving source of the flow (see 
Fig. 4). It should be pointed out, that high velocity gas is also 
present in the vicinity of VLA 1623. This emission is slightly 
offset to the north and aligned with the axis of the outflow. 



3.1.4. Narrow emission feature 

The central emission feature (l<lsr — 2 kms -1 ) is visible in the 
majority of positions in the observed map and diminishes along 
the blue flow axis. The velocity of the peak emission is lowest 
towards VLA 1623 (~ 1 .6 km s ). However, the absorption goes 
deep in this region and no firm conclusions can be drawn regard- 
ing the large-scale velocity field of this component. It is clearly 
affected both by the wing emission (at lower velocities) and the 
deep absorption (at higher velocities). The narrow emission fea- 
ture is strong towards the central position (Fig. 1, spectrum 5) 
and peaks close to the H2CO peak position (Bergman et al. 201 1) 
slightly north from VLA 1623, i.e. at (+17",+76"). 

To the south of the blue-shifted outflow lobe, emission line 
profiles appear single peaked (Fig. 1, spectrum 4). A rough esti- 
mate of the water column density in this region can be obtained 
following the method used by Snell et al. (2000). Assuming ef- 
fectively optically thin gas at a temperature of T = 30 K, we 
estimate the water column density in the (+57", -95") position 
to be 7V(o-H 2 0) = 1.3 X 10 2 %(H 2 ) cirT 2 . 



3.2. PACS map 

The H20(2i2 — lot) line emission also seems to be extended 
along the outflow (see Fig. 2). However, the highest signal-to- 
noise ratio is obtained towards the positions of GSS30, HH313 A 
and HH313 B. In addition, enhanced emission is detected close 
to the position of VLA 1623 (see Sec. 4.1.3). 

3.3. H 2 line survey 

3.3.1. HIFI 

Five H2O transitions were observed with HIFI towards the 
Bl and R2 positions (see Fig. 5 and Table 2). The ortho 
H2O (1 10 — lot) line at 557 GHz (E u /k B — 61 K) is, as expected, 
detected at both positions with a high signal-to-noise ratio. The 
para H 2 0(ln - Ooo) line at 1113 GHz (E u /k B = 53 K) is also 
detected with a reasonable signal-to-noise ratio in both positions 
while the para H 2 (2 n - 2 02 ) line at 752 GHz (EJk B = 137 K) 
has a rather low signal-to-noise ratio at both positions. At both 
positions we observe a difference in flux between the H and 
V polarisations for the latter two lines. The cause of this dis- 
crepancy is at present not understood, but could be due to 
the difference in pointing between the two polarisations and/or 
baseline issues. Consequently, we adopt a somewhat higher 
value on the uncertainty, viz. 40% for the 1113 GHz line and 
50% for the 752 GHz line. The ortho H 2 0(2i 2 - lot) fine at 
1670 GHz (E u /k B = 1 14 K) is not detected towards the R2 po- 
sition but there is a possible detection towards Bl (2<x). The or- 
tho H 2 (3n - 3 03 ) line at 1097 GHz {E u /k B = 249 K) is not de- 
tected in any of the observed spectra. 

3.3.2. PACS 

Four H2O transitions were observed with PACS in two regions 
centred on Bl and R2 (see Fig. 6 and Table 2). Towards the 
Bl position, only the ortho H20(3o3 — 2n) line at 1717 GHz 
(E u /k B = 197 K) is detected with a good signal-to-noise ratio. 
Towards the R2 position, the H20(3o3 - 2 12) line and the para 
H 2 (3 B - 2 02 ) line at 2164 GHz (EJk B = 205 K) are detected 
with a good signal-to-noise. The ortho H2O (22i - 1 10) line at 
2774 GHz (EJkg = 194 K) and the para H 2 O(4 04 - 3 i3 ) line 



6 



P. Bjerkeli et al.: H2O line mapping at high spatial and spectral resolution 



-2.5 km/s 
+7.5 km/s 




-17.5 km/s 
+22.5 km/s 




-32.5 km/s 
+37.5 km/s 





-7.5 km/s 
+12.5 km/s 




-22.5 km/s 
+27.5 km/s 



"Q, 





-37.5 km/s 
+42.5 km/s 



-12.5 km/s 
+17.5 km/s 




-27.5 km/s 
+32.5 km/s 




0.04 -» 1.07 K 




100 -100 

RA Offset (") 



100 -100 

RA Offset (") 



-20 



-10 







10 



20 



Velocity (km s ) 



Fig. 4. Integrated HiOOio - loi) intensities in different velocity intervals where Av = 5 kms~'. The center velocity of the bins is indicated in 
each panel. Contours are in steps of 0.15 K kms~' and the first contour is at 0.1 K kms" 1 for all maps except the two maps at lowest velocity, 
where the first contour is at 0.3 K kms -1 and the contour interval is 0.3 K kms -1 . High velocity gas (l<lsr > 35 kms"') is observed both close 
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represent the positions of the sources presented in Fig. 2. The insets in the lower right corners illustrate each velocity interval for the spectrum 
towards VLA 1623. The lower right panel shows the position-velocity map in the direction connecting the two plus signs. Here, ^lsr = 3.7 km s" 1 
is indicated with a dashed line. 



at 2392 GHz (EJk B = 319 K) are not detected towards Bl and 
R2, however, in the spaxel centred on (+4. 6", +1.5"), emission 
is picked up from the central source. 

4. Discussion 

4.1. Density, temperature, and H 2 abundance 

4.1 .1 . Constraints from H 2 mid-IR emission 

Liseau & Justtanont (2009, and references therein) estimate 
the H2 density in the outflow to ~3 x 10 5 cm 4 , the tem- 



perature to ~10 3 K, and the average column density to 
N(U 2 ) - 3.5 x 10 19 cm" 2 . These authors used ISO-CAM to ob- 
serve the 0-0 S(2) to 0-0 S(7) pure rotational lines of H 2 . We 
can however not exclude the possibility that the water emission 
stems from a region with different excitation conditions. For that 
reason, we have used Spitzer H 2 line data (Neufeld et al. 2009, 
Neufeld, private communication) to generate temperature and 
column density maps of the north-western outflow. From the ob- 
served 0-0 S(0) to 0-0 S(7) lines rotational diagrams are made 
in each pixel of the observed map. These fluxes are corrected 
for extinction, using the Rieke & Lebofsky (1985) curve and the 
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dust opacities are taken from Ossenkopf & Henning (1994, no 
ice mantles). The value of the visual extinction is taken to be 
Ay =21 (Liseau & Justtanont 2009), i.e. equal to the mean value 
presented in that paper. In reality, the extinction varies from 3 to 
50 in the region mapped by ISO-CAM. However, the uncertainty 
attributed to this variation is less than 20% for the derived mean 
temperatures in the north-western part of the flow (see below). 
The rotational diagrams show that the lowest energy transitions 
trace a gas component that is colder than the gas traced by the 
lines observed with ISO-CAM. To derive the temperature of the 
low temperature component, we fit a straight line to the first two 
/7-H2 lines, i.e. the S(0) and S(2) lines. For the hot gas, we fit the 
S(3) to S(7) lines. For the warm gas component (low tempera- 
ture in the rotational diagram) we find a temperature of ~400 K 
while the temperature of the hot gas is estimated at ~1000K. 
The column density of hot H2 (~ 10 19 cm -2 ), obtained from the 
analysis of Spitzer data, is consistent with the findings presented 
by Liseau & Justtanont (2009). The column density of the warm 
component discussed here is, however, significantly higher and 



peaks towards HH 313 A at 7V(H 2 ) - 3 x 10 zu cirT 2 (see Fig. 7). 
The spatial extent of the H2O outflow agrees with the spatial ex- 
tent of warm H2 (see Figs. 2 & 7). 



4.1.2. H 2 excitation 

To interpret the water emission, we have used the non-LTE ra- 
diative transfer code RADEX (van der Tak et al. 2007). A grid of 
models was used, where the temperature was varied from 100 K 
to 2000 K, the density from 1 x 10 4 cirT 3 to 1 x 10 8 ctrT 3 and the 
water column density from 3 x 10 11 cirT 2 to 1 x 10 15 crrT 2 . The 
code was run in LVG (Large Velocity Gradient) mode and the 
collisional rate coefficients were taken from Daniel et al. (201 1). 
Using these collisional rate coefficients, the estimated H2O col- 
umn densities are lower by less than a factor of two, compared 
to when the collisional rate coefficients presented by Faure et al. 
(2007) are used. 
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Fig. 8. Integrated line strengths as a function of density, temperature and 0-H2O column density. The resolution of the PACS (1670 GHz) data 
has been degraded to match the resolution of the HIFI (557 GHz) data. The black, cyan and orange curves show the line ratio as a function of 
increasing o-H 2 column density (from 10 12 to 10 14 cirr 2 ) for three different values on the kinetic temperature of the gas (100, 500, and 1500 K). 
The lowest column density is to the left in this figure. Within each curve the density of the gas varies from 10 5 cirr 3 to 10 7 cirr 3 (from left to right). 
The observed beam averaged values are plotted with coloured squares. Red is towards the position of R2, blue is towards Bl, black is towards 
GSS 30, and green is towards VLA 1623. The dashed lines represent the same regions for smaller source sizes where consecutive squares indicate 
source sizes of 10", and 1". Black and purple dots represent observed ratios towards the region close to GSS 30 and the outflow respectively, i.e. 
where j r m b(557 GHz) du > 2.5 K km s" 1 . The inset in the lower right corner shows the line ratio over the map on a linear scale. Here, the offsets 
are with respect to VLA 1623. 



In Fig. 8, we plot the integrated line intensity ratios for the 
H20(lio - lot) and H20(2i2 - lot ) transitions versus the inte- 
grated line intensity for the H2O (1 10 - lot ) transition, assuming 
a constant line width of 20 km s . For these models, the optical 
depths are low (r « 1) except in the case where the column den- 
sity is 7Y(o-H20) = 10 14 cirr 2 . In this case, the optical depths in 
the two lines are close to unity. When comparing with the obser- 
vations, a common resolution of 32" was used for the 557 GHz 
and 1670 GHz data. We find this to be the optimum resolution 
since, the Bl position is located very close to the central source 
and the resolution of the PACS map needs to be degraded in 
order to obtain clear detections toward the selected outflow po- 
sitions. 

To compare with the observational data we plot the observed 
total integrated intensity ratios in a few selected positions in the 
outflow (see Fig 8). If the beam-filling factor is smaller than 
unity, we assume that the source structure and the beam shape 
are Gaussian (dashed lines in Fig. 8). We note that the total inte- 
grated intensity over the 557 GHz line (including absorption and 
emission) is higher than the corresponding value measured over 
the line wings (excluding absorption and emission) by ~40% in 
the R2 position and by ~20% in the Bl position. 

4.1.3. VLA 1623 and GSS 30 

From Fig. 8 we estimate the beam averaged column den- 
sity towards the positions of VLA 1623 and GSS 30 to 



Af(o-H20) J 1 x 10 13 cm 4 . Assuming that the emission is ex- 
tended over the HIFI beam (coloured squares in Fig. 8), the 
data are consistent with gas of intermediate to high tempera- 
ture (i.e. T > 200 K). We note that the PACS H 2 0(2 12 - l i) 
map does not reveal any extended emission towards the posi- 
tion of GSS 30, i.e. the source seems point-like to Herschel. The 
dashed line represents the intensity ratio for a decreasing beam 
filling factor, down to a source size of 1". In the case where the 
beam-filling factor is small, the data are more consistent with a 
scenario where the column density is higher than 1 x 10 13 cm -2 . 

As mentioned in Sec. 3, enhanced H20(2i2 - lot) emission 
is detected close to the position of VLA 1623. Comparison with 
7 mm interferometer maps from Ward-Thompson et al. (2011) 
shows that the water emission originates in a region that is ex- 
tended in the direction orthogonal to the outflow axis but also 
along the outflow axis (see Fig. 9). The angular extension from 
the position of VLA 1623 is ~8" to the north-east and ~5" to 
the south-west. These numbers correspond to 9 x 10 2 AU and 
6x 10 2 AU respectively. In other words, this feature is larger 
than the measured Full Width Half Maximum (FWHM), 70 AU 
(53 AU at 120 pc), of the inner disk (Pudritz et al. 1996). The 
line intensity ratio of the 1670 GHz and the 557 GHz lines 
allows us to set some limits on the physical parameters. The 
data are only consistent with a scenario where the tempera- 
ture is T > 300 K and the water column density is higher than 
NifyO) = 1 x 10 14 cirr 2 . Assuming that this emission originates 
in an extended ellipsoidal structure, with a size comparable to 



10 



P. Bjerkeli et al.: H2O line mapping at high spatial and spectral resolution 



x 10 ' 




14 12 10 8 6 4 2 -2 -4 -6 -8 -10-12-14 
RA Offset (") 
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Fig. 11. Colour-scale shows the o-H 2 abundance variation with respect 
to H2 in the north-western part of the flow (see Sec. 4. 1 & Fig. 7). The 
positions of R2, Bl, VLA 1623, GSS30, HH313A and HH313B are 
indicated with coloured dots. 



Assuming that the volume density and the kinetic tempera- 
ture are constant in the north-western flow (n(H>) = 10 6 ctrT 3 
and T = 400 K ), we estimate the beam averaged column den- 
sity in all positions of the map from the H20(lio - loi) inte- 
grated line intensities. Comparison with the warm Hb gas ob- 
served with Spitzer thus allows to estimate the water abundance 
in the region (Fig. 11). The water abundance is only presented in 
the region where the estimated temperature is, T = 400+120 K. 
Changing the volume density by an order of magnitude will al- 
ter the abundance by an equal amount in the opposite direction, 
due to the fact that the water emitting gas is sub-thermal and at 
a relatively low density. It should therefore be pointed out that 
this map merely shows the abundance variation, not the abso- 
lute values. The 0-H2O abundance is estimated to be a few times 
10~ 7 with only small variations in the flow. The derived water 
abundance will be discussed further in Sec. 4.3. 



the angular extension, the total water mass has to be at least 
10 23 g (10- 5 M e ). 

4.1.4. R2, B1 and the north-west jet 

Towards the R2 and Bl positions, the beam averaged wa- 
ter column densities are estimated to be in the range, 
7V(H 2 0) = (0.03 - 10) x 10 14 cm" 2 . The temperature is esti- 
mated to be T > 200 K at Bl and T Z 400 K at R2. In Fig. 10 we 
plot the values for NifyO), «(H2) and T, that are in agreement 
with the observed values within the uncertainties. The ortho-to- 
para ratios for H2O are estimated to be 2.0*ij and 2.0*H at R2 
and Bl respectively. The mean values correspond to a spin tem- 
perature of ~20 K. 

The volume densities in the Bl and R2 positions are con- 
sistent with the density inferred from the large [0 1] 63 /urn to 
[Oi]145/mi ratio presented by Liseau & Justtanont (2009). 
Assuming that the FL column density is equal to the value in- 
ferred from the Spitzer H2 data (warm component), we estimate 
the water abundance towards Bl as XifyO) = 10 -7 . 



4.2. Kinematics and energetics 

In this section we estimate the physical parameters of the flow, 
using the combined information obtained from the FL and H2O 
observations. Assuming that the blue and red flows are rea- 
sonably symmetric, these results can be compared with the 
work presented by Andre et al. (1990). The characteristic ve- 
locity of the H2O gas can be estimated from the observation of 
0-H2O at 557 GHz, using the method described in Andre et al. 
(1990, their Eq. 5) and averaging all positions in the mapped 
area. This yields u c h ar = +16kms _1 for the redshifted part 
of the flow (wlsr > 8 kms -1 ) and = -12kms _I for the 
blueshifted part of the flow (ulsr< km s _1 ). The column den- 
sity of warm FL in the redshifted part of the outflow is typically 
N(H 2 ) - 10 20 cur 2 . The extent of this part of the flow is -150", 
while the width of the flow is of the order of 20". Thus the to- 
tal mass of warm H2 in the red-shifted flow is ~1.6 x 10~ 3 M , 
i.e. slightly lower than what was derived from CO observations 
(Andre et al. 1990; Dent et al. 1995). From the flow extent, mass 
and characteristic velocity, other characteristic flow parameters 
can be derived. Here we use an inclination correction equal to the 
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one adopted by Andre et al. (1990), viz. 1 /cos(i), where i = 80° 
is the inclination with respect to the line of sight (Method 1). 
Since the inclination correction is a crucial parameter when de- 
riving the flow parameters, it should be noted that this correction 
may be too large. Furthermore, the dynamical time scale may 
be overstimated by an order of magnitude when using the char- 
acteristic velocity (see Downes & Cabrit 2007). For that reason 
we also estimate the parameters of the H2O and H2 gas using 
the inclination corrections suggested by these authors (Method 
2). In particular we do not apply any inclination correction when 
computing the momentum of the flow. For the kinetic energy, we 
apply a correction of 5 (see Downes & Cabrit 2007, Fig. 4 and 
5) and we use the most recent estimate of the inclination angle, 
i.e. 75°. The dynamical time scale is computed from the maxi- 
mum velocity of the flow, i.e. 30 km s _1 , deprojected by cos(i). A 
comparison between the derived parameters from CO (2-1) and 
H2O is presented in Table 3. Uncertainties in this table are based 
on the assumption that the inclination angle is correct to within 
5°. This yields an uncertainty range, on the characteristic ve- 
locity of the red-shifted gas, of +45 km s~' < if c h ar < +90 km s 
(Method 2). When using Method 1, the larger characteristic ve- 
locity for H2O compared to CO is predominantly due to the shal- 
lower shape of the H2O line profile. The momentum rate is inde- 
pendent of the molecular tracer used, i.e. the force is the same. 
The same applies to the energy and the mechanical luminosity. 

4.3. Observed line profile shapes 

The HIFI 557 GHz map (see Fig. 1) shows a remarkable vari- 
ety in line profile shapes. In most positions of the map, how- 
ever, the line cores have an emission peak that is blue-shifted 
with respect to the absorption feature. This type of line profiles 
can be a sign of infall and is also commonly observed towards 
Class sources (Kristensen et al., A&A in press). In the case 
of VLA 1623, however, infall profiles are observed over a large 
field of view and they could for that reason be due to the veloc- 
ity structure of p Oph A (Bergman et al. 201 1) instead of being 
caused by large scale infall. The absorption feature is centred 
at a velocity, with respect to the LSR, similar to the peak ve- 
locity of the C ls O(3-2) emission (Liseau et al. 2010), i.e. at 
wlsr — 3.7 kms _I . C I8 peaks in the dark core region, close to 
the position where the absorption is strongest in the 557 GHz 
map. It is thus a possible scenario that wide spread H2O emis- 
sion is observed in the line center and the profile suffers from 
most absorption towards the dark core at a velocity, slightly off- 
set. In the velocity interval, 2 to 3 kms , where the 557 GHz 
emission peaks, also wide spread C 18 has been detected (see 
Fig. 4 of Liseau et al. 2010). 

As already discussed in Sec. 3, both blue-shifted and red- 
shifted gas is observed in the north-western region of the map 
due to the presence of other flows. However, the south-eastern 
outflow lobe provides a clean case in the sense that it is free 
from contamination from other flows. From the observed line 
profiles in this part of the flow we can draw some conclusion 
regarding the geometry. The inclination angle of the VLA 1623 
outflow has been estimated at ~75° with respect to the line of 
sight (Davis et al. 1999, from observations of HH313). In the 
extreme scenario where the flow is laminar and only shear is 
present, this would yield maximum velocities of ~ 140 kms~' 
(assuming that the inclination angle of the south-eastern flow is 
comparable to the inclination angle of the north-western flow). 
This is higher than the maximum velocities of ~ 60 - 80 km s 
that have been estimated from proper motion studies (Caratti o 
Garatti et al. 2006; Davis et al. 1999). One could argue however, 




Fig. 12. Cartoon image illustrating a cut through the model described in 
Sec. 4.3. The flow is inclined 75° with respect to the line of sight. In 
this figure, the thickness of each shell has been increased by a factor of 
20 for clarity. The arrow points in the outflow direction and the dashed 
lines indicate the line of sight and the plane of the sky. Note that the 
surrounding cloud is not included in this image. 



that higher velocities than this are present in the jet of the flow. 
However, from Fig. 4, there is no obvious indication that molec- 
ular gas moving at a high velocity is highly collimated. We there- 
fore conclude that transverse gas motions may play an important 
role. Hence, the scenario can be compared to the one presented 
in Bjerkeli et al. (2011) for HH54 (apart from the fact that the 
VLA 1623 outflow lies almost in the plane of the sky). These 
authors conclude that the geometry may be important when it 
comes to the observed line profile shapes. 

In Sec. 4. 1 we put constraints on the H2 density, kinetic tem- 
perature, and 0-H2O column density. To test whether the ob- 
served line profiles can be explained by a curved geometry, or 
not, we use a Monte Carlo code (Hartstein & Liseau 1998) to 
compute the line profiles, where the most recent collisional rate 
coefficients have been used (Daniel et al. 201 1). The model con- 
tains eight spherical shells (see Fig. 12) of a diameter equal 
to the estimated width of the flow, i.e. 4 xl0 16 cm (20" at 
<i=120pc). The thickness of each shell (see Table 4) attains 
the highest value in the center of the flow and goes to zero at 
the edges. The velocity of the gas is varied linearly from 40 
to kms -1 , and the volume density and kinetic temperature is 
set to n(H 2 ) = 1 x 10 6 crrT 3 and T = 400 K respectively. The 
shells are surrounded by a low-density, low-temperature cloud 
with a diameter of 0.3 pc, responsible for the absorption at 
i^lsr = 3.7 kms -1 . 

A synthesized spectrum is presented in Fig. 13, together with 
the observed spectrum towards the B 1 position. In the presented 
model, the H2O abundance is slightly lower than what was in- 
ferred from the observations towards the north-western part of 
the outflow (see Sec. 4.1 and Fig. 11), i.e. X(H 2 0) = 6 x 10~ 8 . 
Other parameters of the model are presented in Table 4. Due to 
the relatively large velocity gradients and the curved geometry, 
the optical depths in the line wing are low (r <K 1), although the 
optical depth at the line center is as high as 1 1. Gas motions per- 
pendicular to the flow axis, therefore, show up as a red-shifted 
emission feature. Further out in the south-eastern flow, however, 
the red-shifted component is less prominent. The spectrum to- 
wards (+76",-38") is thus not easily reconcilable with a sce- 
nario where the inclination angle with respect to the line of sight 
is as high as 75°. In this position, the line profile appears tri- 
angular with no, or little red-shifted emission. The synthesized 
line profile is in better agreement with the observed spectra in 
the north-western part of the flow. It should be noted, however, 
that gas originating from GSS 30 (blue-shifted) is also present in 
this region. 

Previously published H2 (2.12 fim) images reveal a bow- 
shaped structure (Davis et al. 1999) towards the position of 
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Table 3. Comparison between the flow parameters derived from CO (2-1) and H 2 0. The uncertainty on the inclination is assumed to be 5° and 
the velocity of the wind is set to 200 km s~' . 



Parameter' 


Symbol and unit 


Method 1 


a (i = 80°) 


Method 2 b (i = 75°) 


CO (2- 


1) 


H 2 (red) 


H 2 (red) 


Length to width ratio 




14 




7.5 


7.5 


Total mass 


Mtot («g) 


10 x 10~ 3 


t.o X tU 


1 / v 1 n-3 
1.0 X 1U 


Extent 


^rnax (PC) 


0.08 




0.09 


0.09 


Characteristic velocity 


*>ehai (kms~') 


50±?J 




7 -30 


60:11 


Dynamical time scale 


U (yr) 


1-7^ 


x 10 3 


1.0 x 10 3 


7.4 +J; 4 xlO 2 


Momentum 


P(M a kms-') 


6 +0 2 




1 +0 1 

1 -0.05 


0.03" 


Kinetic energy 


E (erg) 


3+ 3 x 


Q 44 


1^4 x!0 44 


2.0 x 10 43 


Momentum rate 


P(M e kms 1 yr 1 ) 


3 4 +4 - 2 

-2.6 


x 10- 4 


1.5 4 _ 5 8 xlO- 4 


3.4i 7 8 xlO- 5 


Mechanical luminosity 


E(L Q ) 


1 4 +34 




1 1 +78 


"•^ -0.04 


Wind mass-loss rate 


M (M a yr~') 


, 7 +2.1 
L -' -1.3 


X 10" 6 


0.7 +oi x 10-" 


1.7 xlO- 7 



Notes. (n) The parameters derived from CO (2-1) were presented in Andre et al. (1990) and have here been scaled to the distance 120 pc. {b) The 
inclination corrections are based on the work presented in Downes & Cabrit (2007). The dynamical time-scale is derived from the maximum 
velocity of the H 2 gas, i.e. ~30km s" 1 . (c) The length to width ratio and total mass are inferred from the Spitzer H 2 observations. 
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Fig. 13. Left panel: Observed 557 GHz spectrum (solid black line) at the B 1 position compared to the spectrum computed with the Monte Carlo 
code (solid red line) where X(H 2 0) = 6 x 10~ 8 . Note that the narrow emission component is not included in the model. Right panel: The observed 
spectrum (black solid line) at the HH313 A position compared to the computed spectrum (red solid line). The model has two components, an 
extended component (dashed line) and a small hot component (dashed-dotted line). Note that blue-shifted emission, originating from GSS 30, is 
visible in the observed line profile. 



Table 4. Parameters used in the Monte Carlo models. 



Parameter and symbol 


Bl 


HH 3 13 A (hot) 


Distance to source, d 


120 pc 


120 pc 


H 2 abundance, X(H 2 0) 


6x 10" 8 


1 x 10- 5 


LSR velocity, i>lsr 


+ 3.7 kms" 1 


+ 3.7 kms" 1 


Velocity profile, v(r) 


40r(0<r< 1) 


40r(0<r< 1) 


Maximum shell thickness 






ratio, A r/R mllx 


5 x 10" 2 


1 x 10" 2 


Source size, 8 SOU me 


20" 


5" 


Microturbulence, v tUT b 


1.0 kms" 1 


1.0 kms" 1 


H 2 density, ;i(H 2 ) 


1 x 10 6 cm" 3 


1 x 10 6 cm" 3 


Kinetic gas temperature, 7\j n 


400 K 


1000 K 



HH 313 A. This is also the position where the water wing emis- 
sion is the strongest and the feature is clearly visible in the 
Spitzer H> data (see Fig. 7). In this region, hot gas is present 
at a temperature of -1000 K, which is why an additional high- 
temperature component is required to reproduce the observed 
emission line profile correctly. The red-shifted part of the emis- 
sion line profile can be explained if an additional shell with a 
size of 5", is added to our previous model. The temperature 
of this emission is taken to be equal to the temperature of the 
hot H2 gas (i.e 1000 K) and the volume density is again set to 
n(H2) = 1 X 10 6 cm -3 , i.e. in close agreement with the models 
presented by Davis et al. (1999). The red-shifted part of the ob- 
served emission line profile is satisfactorily reproduced when the 
water abundance in the hot component is X(H20) = 10~ 5 . 
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We conclude that the derived water abundance in the ex- 
tended flow is less than 1CT 6 . It can, however, be as high as 1CT 5 
locally in the shocked region, HH 3 1 3 A. 

5. Conclusions 

• We have obtained maps of the two ground-state rotational 
transitions of ortho-water, towards the outflow emanating 
from VLA 1623. In addition to this, selected transitions of 
ortho- and para-water have been observed during a line sur- 
vey conducted towards two positions in the outflow. 

• The H2OU10 - lot) and H20(2i2 - lot ) lines are detected 
in the outflow with the HIFI and PACS instruments. In 
addition to this, the H2O (2n - 2rj2), H20(3i2-3o3) and 
H2O (3o3 - 2n) lines were detected towards the positions Bl 
and R2. Towards R2, also the H20(3i3 - 202) line was de- 
tected. 

• Observational data are in agreement with 
LVG models where the column density is 
/V(H 2 0) > 10° cnr 2 towards VLA 1623 and GSS 30 
and /V(H 2 0) (0.03 - 10) x 10 14 cnr 2 towards Bl and R2. 
The temperature of the gas responsible for the H2O emission 
is higher than 200 K in the observed region. 

• Emission at 1670 GHz is detected in a spatially extended 
region close to VLA 1623. These data are only consistent 
with a scenario where the column density is higher than 
JV(H 2 0) 10 14 cm" 2 . 

• The inferred water abundance is lower than what is expected 
from recent shock models. In the north-western outflow it is 
estimated to be lower than 10~ 6 . However, in the HH313 A 
bow-shock on spatial scales of 5", it is estimated to be 10~ 5 . 

• Estimates of the momentum rate, energy and mechanical lu- 
minosity, suggest that the force driving the VLA 1623 out- 
flow is invariant with respect to the molecular tracer used. 

• Line profile shapes in the south-eastern flow suggest that 
transverse gas motions play an important role. The inclina- 
tion with respect to the line of sight is likely smaller than 
75°. 
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